Spectral variability of Ultra Luminous Compact X-ray Sources 

in Nearby Spiral Galaxies 
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ABSTRACT 

Using the X-ray data taken with ASCA, a detailed analysis was made of intensity and spectral 
variations of three ultra-luminous extra-galactic compact X-ray sources (ULXs); IC 342 source 1, 
M81 X-6, and NGC 1313 source B, all exhibiting X-ray luminosity in the range 2 x 10'^^-1.5 x 
10'"' erg s~^. As already reported, IC 342 source 1 showed short-term X-ray intensity variability 
by a factor of 2.0 on a typical time scale of 10 ks. M81 X-6 varied by a factor of 1.6 across 
seven observations spanning 3 years, while NGC 1313 source B varied by a factor of 2.5 between 
two observations conducted in 1993 July and 1995 November. The ASCA spectra of these 
sources, acquired on these occasions, were all described successfully as optically-thick emission 
from standard accretion disks around black holes. This confirms previous ASCA works which 
explained ULXs as mass- accreting massive black- hole binaries. In all three sources, the disk color 
temperature was uncomfortably high at Tin = 1.0-2.0 keV, and was found to vary in proportion 
to the square-root of the source flux. The apparent accretion-disk radius is hence inferred to 
change as inversely proportional to Tin. This suggests a significant effect of advection in the 
accretion disk. However, even taking this effect fully into account, the too high values of Ti„ of 
ULXs cannot be explained. Further invoking the rapid black-hole rotation may give a solution 
to this issue. 

Subject headings: black hole physics — galaxies: X-rays — X-rays: galaxies 



1. Introduction 

One of the most puzzling aspects of X-ray emis- 
sion from normal galaxies is the presence of ex- 
otic, powerful, apparently point-like X-ray sources 
in arm regions of some spiral galaxies (Fabbiano 
1989). We call them "ultra- luminous compact X- 
ray sources (ULXs)", because their luminosities, 
typically 10'^^"^° erg s~^, are one to two orders of 
magnitude higher than those of the most luminous 
Galactic X-ray binaries. Since the Eddington limit 



for an object with mass M is expressed as 



Le = 1.5 X 10 



38 



M 



erg s 



(1) 



where Mq denotes the solar mass, ULXs were 
regarded as mass accreting black-hole binaries 
(BHBs) involving massive (50-100 Mq) stellar- 
mass black holes (BHs), or neutron-star bina- 
ries (NSBs) with highly coUimated X-ray emis- 
sion. However, both explanations remained specu- 
lative, because X-ray spectra of ULXs were poorly 
known. 

The breakthrough has been brought by ASCA 
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(Tanaka ct al. 1994), thanks to its fine spectral 
capability with a modest angular resolution. Fol- 
lowing previous works on individual ULXs (e.g., 
IC 342 source 1; Okada ct al. 1998), Makishima 
et al. (2000) have performed an extensive spec- 
tral investigation of ULXs, and found that their 
spectra arc commonly represented by a so-called 
multi-color disk-blackbody (MCD) model that de- 
scribes emission from an accretion disk around a 
compact star (Mitsuda ct al. 1984; Makishima et 
al. 1986), strongly reinforcing the BH interpreta- 
tion. The MCD model provides two basic quan- 
tities of the accretion disk; the inner-most disk 
temperature, Tin, and the inner- most disk radius, 
Rin ■ Comparing these two quantities of ULXs with 
those of Galactic/Magellanic BHBs, Makishima ct 
al. (2000) pointed out that the obtained values of 
Tin are significantly higher than those of ordinary 
BHBs, whereas the values of _Rin arc not so much 
different. Since i?in might correspond to the last 
stable orbit around a central BH, this fact implies 
that the mass of BHs in ULXs arc relatively low 
(~10-20 Mq), being inconsistent with their high 
luminosities. To cope with this issue, they sug- 
gested that the central BH in a ULX is rapidly 
rotating (i.e., a Kerr BH), hence the disk can get 
closer to the BH and the observed value of Tin 
can become higher. This idea was originally pro- 
posed by Zhang et al. (1997) to explain the similar 
high values of Tin found from the two Galactic jet 
sources, GRO J1655-40 and GRS 1915-1-105. 

Although our understanding of ULXs have thus 
made a great progress, previous works have been 
dealing only with the time-averaged spectra, and 
left their spectral variability unexamined. Such 
a study would bring us key information on the 
source nature. For example, Ginga observations 
of variable Galactic/Magellanic BHBs (Ebisawa 
1991; Ebisawa et al. 1993) revealed that their 
source luminosities vary as oc T^, or _Ri„ re- 
mains constant; this fact implies that an optically- 
thick, geometrically-thin standard accretion disk 
(Shakura & Sunyacv 1973) is actually realized 
around a Schwarzschild BH, where the inner disk 
radius is defined by the last stable circular orbit. 
We also point out that in the "Kerr BHBs" sce- 
nario of ULXs, only a decrease of the inner disk 
radius has been considered. Needless to say, it 
is also necessary to examine how the relativistic 
effects (i.e., gravitational redshift, Doppler shifts. 



and the light bending) affect the observed spectra. 

Accordingly, we performed (re-) analysis of the 
spectra of three luminous ULXs; IC 342 source 1, 
M81 X-6, and NGC 1313 source B. The first has 
the largest flux among the observed ULXs next to 
M33 X-8, and showed strong short-term variability 
(Okada et al. 1998). The other two have relatively 
high fluxes and have been observed several times 
with ASCA, and are hence suitable for the study 
on the spectral variability. We describe the ASCA 
observation of these sources in the next section, 
and examine how their spectra varied in § 3. We 
discuss their source natures in § 4. 

2. Targets and Observations 

IC 342 is a nearby starburst Scd galaxy located 
close to the Galactic plane {b 10°). After TuUy 
(1988); wc assume the distance as 3.9 Mpc, al- 
though large optical extinction caused by the low 
Galactic latitude makes the distance uncertainty 
relatively large. This galaxy has been observed by 
ASCA twice; in 1993 September and 2000 Feburu- 
ary. As already described in Okada et al. (1998), 
source 1 in this galaxy exhibited a significant vari- 
ability during the 1993 observation (see also Fig- 
ure 1). Accordingly, we investigate the change of 
its spectral property in detail in § 3.1. The spec- 
tral difference between the 1993 and 2000 obser- 
vations is described in a separate paper (Kubota 
et al. 2001). 

M81 (NGC 3031) has a Cepheid-based accurate 
distance of 3.6 Mpc from the Hubble Space Telec- 
sopc observations (Frecdman ct al. 1994). Fab- 
biano (1988) observed this galaxy with Einstein 
and detected several sources, including the bright- 
est one, called X-5 (Ishisaki ct al. 1996; lyomoto 
1999) at the nucleus, and the most luminous off- 
nucleus source called X-6 which we study in the 
present paper. Due to the explosion of SN 1993J, 
ASCA has frequently observed the M81 region 
(Kohmura 1994; Kohmura et al. 1994; Uno 1997). 
Among these observations, wc utilize the data ob- 
tained after 1994, in which SN 1993J had faded 
away significantly and its contamination to the X- 
6 spectrum is sufficiently low. Thus we analyzed 
data from seven observations in total, as shown in 
Table 1. Although the first two datasets were al- 
ready investigated by Makishima ct al. (2000), we 
here deal with them as well as other five, in order 
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to study the long-term variability. The same seven 
datasets were also employed by Uno (1997), in his 
study of SN 1993J. 

NGC 1313 is a nearby face-on, late-type Sb 
galaxy at a distance of 4.5 Mpc (Vaucouleurs 
1963) . X-ray observations using Einstein Observm 
ROSAT, and ASCA (Fabbiano & Trinchieri 1987; 
Colbert et al. 1995; Miller et al. 1998; Petre et al. 
1994) showed that its X-ray emission is dominated 
by three extremely luminous point-like sources of 
Lx 10^^ erg s~^ each. One of them is an X-ray 
luminous supernova, namely SN 1978K, and the 
other two are ULXs called source A and source B 
after Petre et al. (1994). The former is located 
close (~45") to the galaxy nucleus, whereas the 
latter is seen at the south-end of the host galaxy. 
This galaxy has been observed by ASCA twice, in 
1993 July and in 1995 November. Makishima et 
al. (2000) studied source B, and showed that the 
spectrum of this source on either occasion can be 
expressed by the MCD model with different values 
of Tin- We here investigate this ULX for both ob- 
servations in § 3.3, concentrating on the spectral 
variability. As for source A, the two observations 
will be compared in another paper (T. Mizuno et 
al. in preparation). 

The observational log of our target sources is 
summarized in Table 1. For M81 X-6, we did not 
use the CIS data (see § 3.2). The CIS data for 
NGC 1313 in the second observation was acquired 
in a bit assignment (PH-X-Y-RT-SP-Time) of 
8 8 8 0-0-7; this is an exceptional mode which 
sacrifices rise time information to improve the time 
resolution (in search for a pulser in SN 1978K), 
so that we did not apply the off-line rise-time cut 
screening on the data. For the other CIS data and 
the SIS data, wo applied a standard data screening 
and tabulated the net exposure in Table 1. 

3. Data Analysis 
3.1. IC 342 source 1 

The time-averaged spectra of IC 342 source 1 
were already analyzed by Okada ct al. (1998) and 
Makishima et al. (2000). They accumulated the 
on-source spectra from a circular region of 3' ra- 
dius on the source center, subtracted the back- 
ground spectra extracted from blank-sky observa- 
tions, and fitted the obtained spectra with sin- 
gle component models. According to their results, 



the power-law model is completely imacceptable, 
the thermal bremsstrahlung (TBS) model provides 
much better value of but is also rejected at 99% 
confidence level, whereas an MCD one gives an ac- 
ceptable fit. 

Although the short-term variability of source 1 
was already sketched by Makishima (1994) and 
Okada et al. (1998), we reproduced its GIS light 
curves in Figure 1. As already reported by Mak- 
ishima (1994) and Okada et al. (1998), the source 
exhibits clear time variability by a factor of 2 in 
a few hours, especially in the hard energy band. 
In order to study this behavior in more detail, we 
divided both the SIS and the GIS data into five 
time regions as defined in Figure lb, and analyzed 
the corresponding time-sorted spectra. To grasp 
the rough information of the spectral variability, 
we first summed the phase 2 and 4 spectra into a 
"low-flux phase" spectra, and summed the phase 
1, 3, and 5 spectra into "high-flux phase" spec- 
tra, as shown in Figure 2. In agreement with the 
light curves, the difference between the two spec- 
tra is more signiflcant at higher energies, indicat- 
ing that the spectrum hardens as the source flux 
increases. We fitted the high- flux phase SIS/GIS 
spectra jointly with the MCD model, and obtained 
the result as given in Table 2; the fit turned out to 
be acceptable {x/v = 129.0/106). We also fitted 
the low-flux phase spectra; the flt is again accept- 
able (x/z^ = 115.6/107). The best flt values of the 
line-of-sight absorption are A^h = 4.7 x 10^^ cm~^ 
for both spectrua, and are consistent with the 
time-averaged one obtained by Makishima et al. 
(2000). We also found that the disk tempera- 
ture differs by a factor of 1.3 between two phases, 
whereas the source flux by a factor of 1.7 

In order to examine changes in the physical con- 
dition of the accretion disk, we c;ak;iilated the val- 
ues of -Kin for the spectra of high-flux and low-flux 
phases. Assuming a face-on geometry, they be- 
came 112 ± 8 km and 149 ± 19 km, respectively, 
where we adopted the ratio of a color temperature 
to an effective temperature of 1.7 after Shimura & 
Takahara (1995), and applied a correction for the 
inner-disk boundary condition after Kubota et al. 
(1998). To our surprise, the radius Rin does not 
appear to be constant but increases as the flux 
decreases; this behavior contradicts that of ordi- 
nary BHBs (Ebisawa 1991; Ebisawa et al. 1993). 
To investigate this inference more quantitatively, 
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we fitted the two sets of spectra simultaneously by 
constraining iVn and i?in to take the same values 
between them (while allowing Tin to vary indepen- 
dently) , and then letting i?in to change separately. 
Then, the fit improved from 'x^ jv = 266.4/216 to 
X^/zv = 245.0/215, indicating that the change of 
i?in is statistically significant (at 99% confidence 
level by an F-tost). 

Although i?in thus changes significantly, it 
might be an artifact due to wrong modeling of 
the spectra. For example, the MCD emission 
from NSBs is usually accompanied by a black- 
body hard component, and BHBs in the soft 
state often exhibit a spectral hard-tail in ad- 
dition to the MCD component. Fluctuation of 
such a hard component, if any, might apparently 
cause the correlated changes of Tin and iiin seen 
from IC 342 source 1. We therefore tried to fit 
the high/low-flux phase spectra simultaneously 
by a common MCD model with the same i?in 
and Tin, adding a hard-component model to the 
high flux phase spectra only. The hard compo- 
nent was modeled by a black-body of temperature 
fixed at 2.0 keV, typical for NSBs (Mitsuda et 
al. 1984), or a power-law of photon index fixed 
at r = 2.5, nominal value for BHBs in the soft 
state (Tanaka & Lewin 1995). Then, the MCD 
plus power-law model coiild not explain the ob- 
served spectra {x^jv = 492.2/216). The MCD 
plus black-body model was also statistically unac- 
ceptable (yX^jv = 312.1/216). Thus, the observed 
spectral variability cannot be due to fluctuations 
of the hard component, but should be attributed 
to the change of the MCD component itself. 

We finally fitted five time-sorted spectra for 
the five time intervals individually by the MCD 
model, with the absorption fixed at the best-fit 
value of high/low-flux phase spectra [N^i = 4.7 x 
10^^ cm^^). We plotted the confidence contours of 
the five spectra on the Tin--Rfn plane in Figure 3, 
where i??^ is proportional to the normalization of 
the MCD model. Then, like we have already found 
in the two phase spectral fitting, i?in gradually de- 
creases as Tin increases, and the relation is approx- 
imately expressed as i?in oc Tj~^. Further exami- 
nation of this result will be discussed in § 4. 



3.2. M81 X-6 

3.2.1. Accumulation of the source and the hack- 
ground spectra 

When studying the spectrum of X-6, we should 
eliminate the contamination from the two nearby 
sources, SN 1993J and M81 X-5, the latter being 
the low- luminosity active-galactic nucleus of M81 
(Ishisaki et al. 1996; lyomoto 1999). This is be- 
cause X-6 is separated only ^1' from SN 1993J, 
and ~3' from X-5, as previously reported by 
Kohmura et al. (1994) and Ishisaki et al. (1996). 

We accumulated the source spectrum of X-6 
from a circular region of 1.'5 radius, to make the 
contamination from X-5 as low as possible. For 
the same reason, we only used the SIS data; the 
poorer spatial resolution of the CIS would increase 
the X-5 contamination. In the obtained spectrum, 
nevertheless, typically ~50% photons still origi- 
nate from X-5 due to its brightness. In order to 
remove this residual contamination, we accumu- 
lated a background spectrum over another region 
having the same size as was used for the on-source 
spectrum. This background region is located op- 
posite to X-6 with respect to X-5, where wc expect 
a similar amount of contamination from X-5. We 
multiplied a constant factor to thus obtained back- 
ground spectrum, and then subtracted it from the 
on-source spectrum. This "scaling factor" was in- 
troduced to take into account the ASCA XRT's 
asymmetric point spread function, and was de- 
termined based on the ray-tracing (Monte-Carlo) 
simulation developed by the XRT team. By sim- 
ulating the X-5 event for each observation (each 
position of X-5 on the focal plane), we can es- 
timate the ratio of photons in the source region 
to the background region. The ratio, or the scal- 
ing factor, is typically 1.5, ranging over 1.2-1.8. 
This procedure also subtracts the non X-ray back- 
ground (NXB) and the Cosmic X-ray background 
(CXB) after multiplying the same scaling factor, 
which is not precise for the NXB and CXB. How- 
ever, we can neglect this effect, since the ratio of 
the NXB -H CXB count rate to that of the con- 
taminated photons from M81 X-5 is estimated to 
be < 10%. 

Thus, we have obtained the spectrum for X-6 
plus SN 1993J. The contamination from SN 1993J 
is difiicult to eliminate because of their short 
separation. Instead, we take it into account 
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in a different way. The SN 1993J spectrum 
was separately estimated by Uno (1997) through 
one-dimensional SIS image analysis, which was 
originally developed by Kohmura (1994). Ac- 
cording to Uno (1997), SN 1993J had faded 
away significantly after one year from the ex- 
plosion, and the spectrum is expressed by a 
power-law of F = 2.5^Qg and the X-ray flux 
/x = 0.5 ± 0.2 X 10 -12' erg s'^ cm'^ in 1994 
April, and of T = 3.0tl ° and fx = 0.2 ± 0.1 x 10 "^^ erg 
in 1994 October, both in the 0.5-8 keV band. We 
calculated their contamination to the X-6 spec- 
trum based on the ray-tracing simulation, and 
took the results into account as a fixed power-law 
component when fitting the X-6 spectrum, as in- 
dicated in Figure 4. The contribution of SN 1993J 
to the 0.5-10 keV flux turned out to be small; 
only ~ 10% and ~ 3% in 1994 April and Octo- 
ber, respectively. For the other data which were 
obtained after 1995 April, we neglected the contri- 
bution from SN 1993 J, since the supernova further 
faded away. 

3.2.2. Spectra of individual observations 

We fitted the obtained seven spectra separately 
with the typical single component models; the 
power-law, the TBS, and the MCD model. Then, 
the MCD model turned out to be always accept- 
able at 90% confidence level, while the TBS model 
provided a worse fit except for the 1995 Octo- 
ber and 1996 October spectra (for the worst case, 
X^/v = 69.8/54 in 1996 April). The power-law fit 
gave the largest values of (1996 April data gave 
the worst fit of x^/t^ = 80.3/54). Therefore we 
conclude that the spectra of M81 X-6, like those 
of IC 342 source 1, can be represented well by the 
single MCD model, and tabulate the MCD-fit re- 
sults in Table 2. Our results for the 1994 observa- 
tion are consistent with those previously reported 
by Makishima et al. (2000). 

We plotted the obtained values of the bolomct- 
ric flux of the accretion disk f^^^ and Tin in Fig- 
ure 5. It indicates that the seven datasets can 
be grossly grouped into two representative states; 
a low-temperature state observed in 1995 Octo- 
ber and 1996 April, and a high-temperature state 
observed on the other occasions. The low/high- 
temperature states correspond to the low/high- 
flux states, so that M81 X-6 shows similar spectral 
variability to IC 342 source 1. We also present the 



line of .fbo^'^ i.e., the locus of Ri^ remain- 

ing constant. Although the statistical errors are 
somewhat large, a slight deviation from this line 
can be inferred. 

We hence made the Ti^-Rf^ diagram in Fig- 
ure 6a, as was done on IC 342 source 1. Thus, 
we again, find a hint of anti-correlation between 
Tin and Rf^, like in the case of IC 342 source 1, 
but relatively large errors hampers definite state- 
s m§§t. We therefore grouped the seven datasets 
into two spectra, i.e., the high- temperature state 
spectrum and the low-temperature state one. We 
fitted these two spectra by the same single com- 
ponent models as were used for each observation, 
and, again, the MCD model turned out to be the 
best representation of the data. We summarized 
the best fit parameters of the MCD model in Ta- 
ble 2, and plotted the Tin Rfn diagram in Fig- 
ure 6b, by fixing the absorption at the average 
value of the two states {Nn = 1.8 x 10^^ cm-^). 
Thus, i?in increases marginally as the flux de- 
creases, although not so noticeable as in the case 
of IC 342 source 1. In fact, when we fitted the two 
subgrouped spectra simultaneously by constrain- 
ing iVfj and -Rin to be common, wc obtained an 
acceptable fit with x2/i^=140. 4/126, while letting 
i?in to be free does not improve the fit signifi- 
cantly (x2/i^=138. 0/125). Therefore i?in is consis- 
tent with being the same between the two spectra, 
although a weak anti-correlation between Tin and 
Rin may be present. 

3.3. NGC 1313 source B 

As previously reported by Makishima et al. 

(2000), this source showed time variability by a 
factor of ~ 2.5 between the 1993 and the 1995 ob- 
servations. They analyzed both spectra with sin- 
gle component models (the power-law, the TBS, 
and the MCD model), and found that the MCD 
one gives the best description of the data. More- 
over, they pointed out that the disk temperature 
is positively correlated with the source fiux, and 
the value of J?in seemed to increase for two years 
whereas the source flux decreased. This behav- 
ior is therefore similar to that found with the two 
sources we have described so far, and we investi- 
gate the variability of this ULX in further detail. 

We here utilized the same data as used by Mak- 
ishima et al. (2000), and fitted them with the 
MCD model. We plot the confidence contours 
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on the Tin-Rf^ plane in Figure 7, where we fixed 
the absorption at iVn = 7 x 10^° cm"^ (the av- 
erage value between the two observations). Like 
the two variable sources so far studied, a slight 
increase in i?in can be seen as the flux decreases. 
The change of Rin is statistically significant (at 
90% confidence); we obtained an improved fit 
(X^/i^=214.9/207) when fitting the two observa- 
tions simultaneously allowing Rjn to take separate 
values, compared with that (x^/i^=222. 9/208) ob- 
tained when constraining R^ to be common. 

Since the source positions on the focal plane 
are largely separated (~ 6 ) between the two ob- 
servations, the obtained value of Rin might be af- 
fected by the response uncertainty. According to 
an extensive investigation utilizing tlic^ Crab Neb- 
ula by Fukazawa et al. (1997), the calculated nom- 
inal source flux may increase artificially as the 
off-axis angle increases. As for our two observa- 
tions, source B was observed at an off-axis angle 
of 10 and ~ 5 in 1993 and 1995, respectively. 
Therefore, if this residual response uncertainty is 
present, the true 1993 flux should be lower than 
the estimated value. This would further enhance 
the increase in i?in for two years. We therefore 
conclude that the change of i?i„ seen for NGC 1313 
source B is a real effect rather than an artifact. 

4. Discussion 

4.1. Two distinctive properties of ULXs 

So far, we have been investigating the time 
variability of ULXs spectra, and found that they 
can always be represented by the MCD model, 
event when the source flux varies significantly. 
We also found that the observed flux is posi- 
tively correlated with Tin. If ^^in remains con- 
stant as the source flux varies like ordinary Galac- 
tic/Magellanic BHBs (Ebisawa 1991; Ebisawa et 
al. 1993; Kubota 2001), we will have the relation 
of /boi'^ oc T^, whereas our ULXs show somewhat 
different relation (e.g.. Figure 5). As summarized 
in Figure 8, they seem to obey a single common 
scaling of i?in oc Tj~^. 

The observed values of Tin of our ULXs, rang- 
ing 1.0-2.0 keV, are also different from those of 
normal BHBs (typically 0.5-1.2 kcV; e.g. Tanaka 
and Lewin 1995). Such a contradiction has been 
known to be a common feature of ULXs (Okada 
et al. 1998; Mizuno et al. 1999; Makishima et al. 



2000). As previously described by Makishima et 
al. (2000), in the MCD approximation for the stan- 
dard accretion disk around a Schwartzshild BH of 
mass M, Tin is expressed as 

--■K^)""(i^)^-(i4)"""- 

(2) 

where r] denotes the disk bolometric luminosity 
normalized to the Eddington luminosity, and we 
have been assuming ^ = 0.412 and k = 1.7. 
Thus, a heavier BH should show a lower value of 
Tin, as expressed by Tn oc M~^/^. Nevertheless, 
the ULXs that have high luminosities and hence 
high BH mass actually exhibit higher disk tem- 
peratures. The upper-limit on Tn inferred from 
equation 2 is 0.68 keV for the phase 1 spectra of 
IC342 source 1, 1.0 keV for the high- temperature 
phase spectra of M81 X-6, and 0.89 keV for the 
1993 spectra of NGC 1313 source B even assuming 
a face-on geometry; these predictions contradicts 
the observed values by a factor of 2.9, 1.6, and 
1.7, respectively. A higher disk inclination angle 
makes the situation even worse. Even if we exclude 
IC 342 source 1, for which the distance to the host 
galaxy is somewhat uncertain, the contradiction 
by a factor of ~ \/3 remains to be solved. 

Thus, our research with ASCA has provided 
ambivalent results on the interpretation of ULXs 
in terms of optically-thick standard accretion disks 
around massive stellar-mass BHs. On one hand, 
we have confirmed that the MCD model remains 
a good representation of the ULX spectra even 
when the source flux varied considerably. On the 
other hand, the innermost disk radius apparently 
changes in an anti-correlation with the disk tem- 
perature, which makes another discrepancy in ad- 
dition to the previously pointed-out problem of 
too-high values of Tin (Makishima et al. 2000). In 
order to solve these two problems, we may intro- 
duce some modification to the standard disk pic- 
ture. 

4.2. Slim disk scenario 

Although the standard disk model successfully 
explains the spectra from BHBs in the soft state, a 
significant progress has been achieved on the the- 
ory of accretion disks. When accretion rate M 
is high and the source luminosity approaches Te, 
the standard disk is predicted to change into a so- 
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called optically-thick, advection-dominatcd accre- 
tion flow (ADAF). This "optically-thick ADAF" 
model (Abramowicz et al. 1988; Szuskiewicz et 
al. 1996; Watarai ct al. 2000) is also named a 
"slim accretion disk model" , since it is moderately 
geometrically-thick. In ULXs, because of their 
high luminosities, the accretion flow configuration 
is expected to change from the standard disk to 
the slim disk. 

One of the most characteristic features of 
the slim disk, revealed by Watarai et al. (2000) 
through their numerical calculations, is that the 
X-rays are radiated not only from the regions out- 
side the last stable orbit, but also from the regions 
inside it, since an abrupt change of the radial in- 
fall velocity docs no longer occur. They fitted the 
numerically calculated slim disk spectrum by the 
MOD model (with i?;^ and Ti„ being the model 
parameters), and found that as M increases, J?i„ 
decreases as Ri^ oc T-^^^. Thus, our finding of 
the change in i?in can be explained naturally by 
presuming that the slim disk is realized in ULXs. 
This suggests the presence of a slim disk in the 
three ULXs studied here. 

Because the slim disk has a smaller i?in and 
a higher Tin than a standard disk in the same 
condition, it could also explain the "too-high Tin" 
problem described in § 4.1. For the present three 
ULXs, the slim-disk scaling (the bolometric lumi- 
nosity Tboi oc Rf^T^ oc Tj^) holds over a typical 
range from the observed highest luminosity Z/max 
(which we tentatively identify with Le) down to 
~ Lniax/2. While the source luminosity varies by 
a factor of 2, we expect the disk temperature to 
change by a factor of 2^/^ for the standard disk, 
or 2^/^ for the slim disk. In this way, the slim 
disk scenario can relax the too- high Ti„ problem at 
least by a factor of ^y/t ^ 1-2- However, in order 
to fully resolve the factor discrepancy seen for 
our three ULXs, the source must make a transition 
from the standard-disk to slim-disk regimes at still 

lower luminosities, e.g., iboi = -ft-^'E ^ 0.1 Le- 

v3 

This contradicts the observed results of Galac- 
tic/Magellanic BHBs, where the standard accre- 
tion disk picture has been confirmed to be valid 
at least up to ~ | Le in case of GS 2000-1-25 and 

LMC X-3 (Ebisawa 1991; Ebisawa ct al. 1993; Mi- 
neshige ct al. 1994; Makishima et al. 2000). 

Through the present study, we suggest that the 



ULXs are in the slim-disk condition because of the 
i?in oc T;~^ property, and that this partially ex- 
plains the too-high values of Tin- However, we at 
the same time presume that the issue cannot be 
fully solved even assuming the slim disk hypothe- 
sis. 

4.3. Spinning BH scenario 

Before the present paper, some authors (Mizuno 
et al. 1999; Makishima et al. 2000) tried to explain 
the "too-high Tin" of ULXs by assuming that the 
central BH is rapidly rotating; this "Kerr-BHB 
hypothesis" has first been proposed by Zhang et 
al. (1997) to explain the observed high tempera- 
tures of Galactic jet sources GRS 1915-1-105 and 
GRO J 1655-40. In this subsection, we examine 
whether the remaining inconsistency seen for our 
ULXs can be solved by considering Kerr BHs or 
not, taking into account the relativistic effects 
that have not been considered by Makishima et 
al. (2000). Hereafter, we express the BH angular 
momentum J in a dimensionless manner, i.e., by a 
spin parameter a* = of^iJ- This parameter takes 
values between —1 to 1; a* = 1 means the extreme 
Kerr hole for a prograde disk (i.e., rotating in the 
same direction as the BH), a* = — 1 also repre- 
sents the extreme Kerr hole but for a retrograde 
disk, and a* = corresponds, of course, to the 
Schwarzschild BH. 

The most immediate effect of the BH spin is 
that it affects the radius of the last stable orbit, 
Riast- While i?iast is 3i?s for a Schwartzshild BH, 
it reduces down to |i?s for a prograde disk around 
an extreme Kerr hole of a* = 1 (Berdeen et al. 
1972). A smaller i?in leads to a higher Tin, sug- 
gesting that the BH spin can explain the problem 
with ULXs. 

Neglecting relativistic effects for the moment, 
we can perform simple quantitative estimates. For 
an extremely Kerr hole, i?in decreases by a factor 
of 6, or Tin at the luminosity maximum increases 
by a factor of y/6. Even taking a somewhat less 
extreme case of a* = 0.95, we expect i?in to de- 
crease by a factor of 3 (or, Tn increases by a factor 
of \/3) compared to the case of a Schwarzschild 
BH. This is apparently sufficient to explain the 
observed high temperatures of ULXs. 

Of course, the X-ray spectra emergent from an 
accretion disk is subject to several relativistic cor- 
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rections, due to gravitational rcdshift, transverse 
and longitudinal Doppler shifts, and gravitational 
focusing. To a distant observer, both the observed 
color temperature and flux will deviate from the 
local values, depending on the inclination angle 
i. These effects have been examined through nu- 
merical calculations by several authors, including 
Cunningham (1975), Asaoka (1989), and Zhang 
et al. (1997). Employing the numerical calcula- 
tion by Cunningham (1975), Zhang et al. (1997) 
represent the relativistic effects by two correction 
factors, which are to be applied on the MCD mod- 
eling of the accretion disk around a Kerr BH. One 
correction factor is the change of the color temper- 
ature denoted kgr, due to the gravitational red- 
shift and Doppler red/blue-shifts. The other is the 
change of the flux denoted g, due to the viewing 
geometry and the gravitational focusing. The ob- 
served color temperature scales as cx kgr, whereas 
the observed flux scales as fx oc ^ instead of 
fx oc in the Newtonian case, where i denotes 
the disk inclination angle. The latter correction 
affects the inferred upper-limit of the disk tem- 
perature as (X i^o/^ oc g^^^. Further taking into 
account the correction due to the decrease of i?in, 
which increases the maximum disk temperature as 

— 1/2 

OC , we calculated the combined correction 

factor, KGjig^^'^R^^^ , for a nearly-extreme Kerr 
holes as shown in Figure 9. 

Thus, if we assume a Kerr BH of a*=0.998, and 
view the disk from an inclination angle of i > 65°, 
the correction factor exceeds -\/2. When this is 
combined with another factor of ^ 1.2 due to the 
slim-disk property (§ 4.2), we can explain the dis- 
crepancy by a factor of -\/3 seen for ULXs. Since 
i ^ 60° is what is expected on average when the 
disks are randomly oriented, the chance probabil- 
ity of finding such objects is reasonably high. We 
thus conclude that the mystery of ULXs may be 
solved in terms of the slim accretion disk and the 
BH rotation. 

Of course, there remain many aspects of the 
ULXs phenomena to be examined. Observation- 
ally, an optical identification of the counter-part is 
important; it can help us to understand how ULXs 
are formed, and confirm whether they are really 
massive, rapidly-rotating BHs. We also need to 
obtain spectra with much higher quality over a 
wider bandpass, in order to compare them with 



the theoretical predictions in detail. Theoretically, 
numerical models of the X-ray spectrum emergent 
from an accretion disk around a Kerr BH, taking 
into account both the relativistic and advection 
effects, are required. Another big theoretical issue 
is how to make massive (^ 100 Mq) and rapidly 
rotating (a* > 0.95) stellar-mass BHs. 
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Time(3) Time(3) 



Fig. 1.— The ASGA GIS2 + GIS3 light curve 

of IC 342 source 1 inckiding background. Panel 
(a) represents the light curve in the energy range 
of 0.7-2.0 keV, whereas panel (b) in 2.0-10 keV. 
The background count rate for the low and high 
energy range is about 0.004 c and 0.006 c s~^, 
respectively. 



Confidence Contours 
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Tin (keV) 



Fig. 3. — The 68% confidence contours of the 
short-term variations of IC342 source 1, in terms 
of Tin and the MCD normalization (proportional 
to The solid line represents the relation of 

Tin oc -Rin^. The horizontal axis is logalithmic. 
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Fig. 2. — The high and the low flux phase spectra 

of IC 342 source 1, fitted with the MCD model. 
The histograms show the best fit model and the 
crosses represent the observed spectra. The lower 
panel indicates the fit residuals. Although the fit 
is simultaneous to the SIS and the GIS spectra, 
here only the GIS spectra are shown for clarity. 



d 
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channel energy (keV) 

Fig. 4.— The SIS spectra of M81 X-6 obtained 
in 1994 April, fitted with an MCD model. Dotted 
lines indicate individual contributions of SN 1993 J 

(lower one) and X-6 (higher one), with the former 
represented by a fixed power-law component (see 
text). 
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Table 1 

asca observational log of the host galaxies. 



Galaxy 


Date'' 


Exposure (ks)'' 


SIS Mode'^ 


SIS Clock^ 




yymmdd 


SIS 


GIS 






IC 342 


93 09 19 


35.8 


38.4 


F/B 


1111/1111 


NGC 3031'' 


94 04 01 


31.6 




F/F 


1000/1000 


(M81) 


94 10 21 


37.4 




F/F 


0010/0010 




95 04 01 


17.6 




F/F 


1000/1000 




95 10 24 


34.5 




F/F 


0010/0010 




96 04 16 


43.2 




F/F 


1000/1000 




96 10 27 


27.8 




F/F 


0010/0010 




97 05 08 


41.0 




F/F 


1000/1000 


NGC 1313 


93 07 12 


21.2 


27.7 


F/B 


1111/1111 




95 11 29 


33.2 


32.9 


F/F 


0100/0100 



^ Observation start date. 

^ An average of the two SIS sensors or the two GIS sensors, after 
the data screening is apphed. 

" Data acquisition mode for Bit-High/Bit-Medium. 

Clocking mode of SO for Bit-High/Bit-Mcdium. 
^ The GIS data are not utilized in the present study (see § 3.2). 
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Table 2 

Time-resolved ULX spectra fitted by the MCD model. Errors are calculated for 90% 
confidence for one interesting parameter (ax^ = 2.7). 



Source 


(1022 cm-2) 


T- 

in 

(keV) 






IC 342 source 1 










high-flux 
low-flux 


0.47 ±0.04 
0.47 ±0.07 


1.87±0.07 
1.42±0.09 


15.5 
8.95 


129.0/106 
115.6/107 


M81 X-6 










1994 April 

1994 October 

1995 April 

1995 October 

1996 April 

1996 October 

1997 May. 
high-temp phase 
low-temp phase 


0.19 ±0.06 
0.23 ±0.06 

0.08 ±0.08 

33+018 
^•"^■^-0.14 

0.08 ±0.07 
0.10 ±0.07 
0.16 ±0.03 
0.20 ±0.08 


1.63±0.19 
1.51±0.14 

1.73±0.27 

1 ..4+0.2.3 

1 9. ,+0.18 
-^•^"^-0.15 

1 63+0-23 
-"^•"•J-o.io 

1.59±0.09 
1.29±0.13 


3.61 
4.05 
3.71 
2.48 
2.78 
4.14 
3.80 
3.92 
2.63 


43.8/42 
34.9/47 
29.6/23 
17.7/21 
61.6/54 
30.2/31 
44.5/40 
68.6/65 
89.8/74 


^ Bolometric 


flux of the 


accretion 


disk in 


units of 



10 ^2 erg s ^ cm assuming a face-on geometry. 
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bolometric disk flux (10 erg s cm 



Fig. 5. — Relation between the disk bolometric 
flux and Tin for M81 X-6. 






1,25 1.50 1.75 

Tin (keV) 



Fig. 8. — Relation between Rm and Tin for the 
three ULXs. The error bars represent the 90% 
confidcmcc errors of the spectral fit, and the solid 
line indicates the relation of i?in oc T-TT^. 



Fig. 6. The same as Figure 3, but for M81 X-6 
instead of IC 342 source 1. (a) The contours of 
individual observations, where only the 68% con- 
fidence levels are shown, (b) The contours of the 
high/low-temperature state spectra obtained by 
grouping the individual datasets. The 68%, 90%, 
and 99% confidence levels are presented. 

Confidence Contours 
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Tin (keV) 

Fig. 7. — The same as Figure 3 and 6, but for the 
NGC 1313 source B, where 68%, 90%, and 99% 
confidence levels are shown. 
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Fig. 9. — Expected increment factor of the disk 

temperature due to a BH rotation, relative to the 
Schwartzschild BH of the same mass and the same 
mass accretion rate. It was calculated based on the 
correction factor of Zhang ct al (1997), for the ac- 
cretion disk around a Kerr hole with a* = 0.9981. 
This value corresponds to the case where a BH 
equilibrates with the disk accretion (Throne 1974). 
The disk is assumed to be in the regime of stan- 
dard accretion disk, with no advection eflFects. The 
dashed horizontal line indicates the goal necessary 
to explain the high Tin of ULXs. 
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